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ABSTRACT 

The standard methods of projecting energy spectra from Hartree- 
Fock intrinsic wave functions become very time-consuming as the 
model space is enlarged. One of these methods, however, can be 
revised so as to minimize storage and handling requirements. A 
careful description of the technique is presented for the general case, 
together with sections concerning simplifications resulting from cer- 
tain nuclear symmetries, and suggestions related to efficiency and 
accuracy in the actual calculation. 

INTRODUCTION 

The abundance of experimental information supporting the shell 
model of the nucleus has made it our most valuable tool for investi- 
gating nuclear structure. Shell model calculations become quite com- 
plex for deformed nuclei, however, because of the large number of 
configurations which must be employed. As a result of this, the 
Hartree-Fock model has in recent years become a very popular sub- 
stitute for the shell model. In some fashion which is not yet entirely 
clear, the extra degrees of freedom afforded by the intrinsic Hartree- 
Fock wave function contrive to produce extrinsic, or physical, states 
which are in many ways equivalent or superior to shell model states. 

One problem associated with the Hartree-Fock model is that 
physical states must be extracted from the intrinsic wave function by - 
some form of angular momentum projection. When one is dealing only 
with expectation values of single- particle operators the problem is not 
too difficult, and has recently been treated in a systematic manner. 
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To find the energy spectrum of a nucleus, however, one needs the ex- 
pectation value of a two-particle operator - and the problem is then 
vastly more complicated. 

Four different methods of energy projection are discussed by 
Ripka.^ Of these, the first two are hopelessly impractical for nuclei 
with more than a few active nucleons. The third method makes use of 
a complete set of intermediate two- particle two-hole states to evaluate 
the required matrix elements; simplifications which arise for even- 
even N = Z nuclei are described in some detail by Tewari and 
Grillot. ^ The basic difficulty with this method is that the sum over 
intermediate states is really a sum over unoccupied orbits, and thus, 
when the model space is enlarged to reduce truncation error, the 
number of terms becomes inordinately large. 

By making use of closure one may transform the sum over un- 
occupied orbits into a sum over the occupied orbits, leading to the 
fourth method mentioned briefly by Ripka. For large model spaces 
this technique should be the most efficient, if properly applied. A 
version employed by Gunye and Warke^ has the disadvantage that it 
involves two-body matrix elements depending on the magnetic quantum 
numbers; the total number of matrix elements which must be handled 
is thus very large. The version presented here constructs coupled 
orbital matrix elements from the coupled two-body matrix elements, 
thus minimizing storage and handling requirements. A complete 
treatment is given for the most general" case, followed by sections de- 
scribing the simplifications resulting from certain kinds of symmetry. 
A concluding section suggests some means of improving efficiency 
and accuracy in the actual calculation. 

THEORY 

In the axially-symmetric Hartree-Fock model for deformed 
nuclei, the intrinsic state of a nucleus is represented by a Slater de- 
terminant <3> which is an eigenstate of J , 

z 

J $ = KK$. 
z 
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but not of J . Physical states of the nucleus must therefore be ob- 
tained by means of a proj ection operator, 


tf, 


«JM 



2J + 1 
2 



d MK (e )* (e)sln 


rsj T 

where $(6) = exp(-i0Jy)<3>. By construction, the operator picks 

out that part of a function with angular momentum quantum numbers J 
and K, and changes it so that its J -eignevalue is RM. 

An important nuclear property is its energy spectrum, which may 
be determined by evaluating the matrix element 
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or more generally, 





where is another (possibly identical) Slater determinant. 

As is often the case, a straightforward approach proves to be the 
most efficient in the final analysis. To begin with we have 



2J + 1 
2 


r d£ tK (0)H(0)sin eae, 

Jo 


H(e) = <$’|h|£> 



= t (e) + v(e). 


where 
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Now the function <f> may be expanded in either of the forms 


*(12 . . .A) = A" 1/2 ^ |N a > * N (23 . . .A) 


= [2A(A - 1)]~ 1/2 J 


N a N b 


l N a N b> *N N h (34 • ' - A) ’ 

a b 


where |N a > is a single-nucleon orbital function with the properties 

J z |Na>=Hm a | N a>, 


<N b l N a> - « N N . 

a b 


and |n N^) is a normalized, antisymmetric orbital pair function: 


i N a N by = ^(l N a>l N b> ‘ I N b> I N a>) 


Using these expansions we obtain 


T(0) = 2, < N i T N a >\*N'- * 


N’N 
a a 


n ' n 

a a 


v(e) = 7 I 2 


4 N a b 1 1 a b' \ N’N/ 1 N N, 

4 N N. N’N’ ' a b a bj 

a b a b 


<N1N’ V N N,) *’ 1Mf $ XT M , 


where the rotated orbitals N 0 ) are defined by 

a 


-iej 

|N a > = e y|N >. 
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Let us now indicate explicitly the dependence of the various matrix 
elements on the rotation angle 6 . We begin with the identity 

K> ■ I P M M l N a>’ 

JM 

from which it follows that 

JMM’ 


A similar relation holds for the orbital- pair state |N and in both 

cases the inner sum over M reduces to one term because the operand 
states are eigenfunctions of J z> The orbital kinetic energy matrix 
element thus becomes 



and, with a similar abbreviation, 

< N a N b |v|Sf a S b > =2 d M ib M ab («< N i N bl V jl N a N b>> 

J 

where M, = m a + m^ and = m^ + m^. We refer to (N^| Tj |N a ) 
and < I Vj | N a N^> as coupled orbital matrix elements. 

The angular dependence of the various overlap integrals cannot be 
displayed so explicitly. It turns out, however, that they are all pro- 
portional to the overlap integral 


A(0)= <$'|$(0)> 
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which is a very sharply peaked function of its argument. This de- 
pendence may be factored out by defining the one- and two-nucleon 
density matrices 




N i N i> 


$ 


N a N b 


\= Me)<N^|p (1) |N a > 

\= A(e){N^|p ( 2 ) |N a N b >, 


in terms of which 


H(e) = A(e)h(e), 
He) = t(e) + v(e), 


*<«> = I I Cm («XN i |T J |N > ><N;|p< 1 ) |M a >, 


N a N a J 


X I 2 d M kb M ab ( e )< N ; N bl V jlW< N a N blP <2) l N a N b >. 


n n, mm j 

a b a b 


The matrices and are angle- dependent, but weakly so, and 
satisfy the relationship 


< N a N bl p(2> l N a N b> = < N ;l p<1) l N a ><Ni,|p (1) |N b > 

- <N;|p ( 1 ) |N b > <N b |p^^|N a >. 

The one-nucleon density matrix p^ is numerically equal to the trans- 
pose inverse of the matrix G whose elements are 


<n: 


, , / J 

-iflJ’ 

|G|N a >= N ; : 

e y 


N 


a. 


I 




^)< N ;l G jl N a>’ 
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and one may further show that A (0) = det {<N^| G|N a )} . 

Next we turn to the question of isotopic spin. The single-nucleon 
orbitals will be assumed to be eigenstates of the z- component of the iso- 
topic spin operator, so that 


The kinetic energy integrand t(0) will therefore divide into two parts, 
tp(0) and t n (0). The first consists of all those terms in t(0) for which 
both N; and N refer to proton orbits; the second consists of the re- 
maining terms, for which both and N a refer to neutron orbits. 
There are no cross terms. 

The potential energy integrand v(0) divides into three parts, v ( 6 ), 
v n (0), and v (0). The first consists of those terms in v(-9) for which 
all four labels refer to proton orbits; the second consists of those terms 
for which all four labels refer to neutron orbits. The third part Vp n (0) 
consists of the remaining terms in v(0), which involve only orbital pairs 
in which one nucleon is a proton and the other is a neutron. 


Using symmetry arguments one may restrict the summation for 
v pn (0) so that and N a refer to proton orbits and and to 
neutron orbits, provided the result is multiplied by four. Similar argu- 


ments allow one to restrict the summation for v (0) and v (0) so that 

p n ' 

n b < n a and n^ < n a , again provided that the result is multiplied by four . 

It is convenient to introduce an isospin projection operator , and 
use the identity 
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for this allows part of the dependence of the coupled orbital matrix ele- 
ments on the labels t , r u to be isolated: 

a * b 

<N'N> | Vj |N a N b > = 2 A T<Vb)k N bl V jl N a N b) T - 

T=0 ' ' T 

The isospin- coupled orbital pair state |N a N^j is defined by 

l N a N b> T “ ^ (jv l“b> + (- 1 > T I V l n a>)'*T’ 

X T being a singlet or triplet isospin function, and the coefficient 

^T^ T a T b^ has the val - ue ^T1 for the inte S rands v p( 0 ) and v n (0) ^ the 
value 1/2 for the integrand Vp n (0). 

Finally we consider the question of parity. It is generally assumed 

that each single- nucleon orbital |N ) is an eigenfunction of the parity 

a 

operator with eigenvalue 7 r; this permits further restrictions on the 

a 

summations. In particular, the orbital pairs (N a ,N b ) needed to evaluate 
v{6) can be classified as shown in Table I: 


TABLE I. - (N a ,N b ) ORBITAL 
PAIR ORDERING 


Integrand 

Pair restrictions 

IsoSpin 

Parity 

Other 

ye> 

(p,P) 

(+,+) 
(+» “) 

% <n a 

n b <n a 

v«> 

i 

1 

(n,n) 

(+> +) 
(+,-) 

n b <n a 

n b <n a 

! 

i 

v (e> 

1 

1 

; 

(p,n) 

(+» +) 

’(+,-) 

(-,+) 
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BASIS FUNCTIONS 

The most general solution of the Hartree-Fock equations would, of 
course, yield the exact radial and angular dependence of all the orbital 
functions. Usually the complexities of this task are avoided by assuming 
that each orbital function is accurately described by an expansion 


l N a> =2 Wl km a>l T a> 
k 

in terms of a limited number of orthonormal (single-particle) basis 
functions I km ) which are eigenfunctions of j with eigenvalue Rm , 
eigenfunctions of the parity operator with eigenvalue ti , and are other- 

w a 

wise labeled by the index k. When the basis functions are also individ- 

2 

ually eigenfunctions of j they are said to form a spherical basis; this 
is the only kind considered here. We assume that the simplified 
Hartree-Fock equations have been solved and that the coefficients C^N^ ) 
are known. 

The elements of the matrix G J5 from which p^\ p^\ and A are 
eventually to be constructed, may now be written 



The coefficients I> k (N^) are the analogues of C k (N a ) for the final- state 
orbitals |N’ > , which we assume are also described by the same set of 

ci f 

basis functions \Jkm a )}. The kinetic energy matrix elements are simi- 
lar in structure, having the form 
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(kgllTllk^) being the reduced matrix element (which for the rotationally 
invariant kinetic energy operator is numerically equal to (kginl T^m) 
for any m). 

To evaluate the potential energy matrix elements, we introduce 
basis -pair states coupled to good J and T, 

|k l k 2>jT )= 2 <h m l’3 2 m 2l JM >K m l k 2 m 2> T 
m l m 2 


in terms of which 



The exchange symmetries of the coupled basis-pair states and the 
Cl ebsch- Gordon coefficients permit this to be written in the abbreviated 
form 


Vb) T = I 




< M ab> 
JT 


TT 1 

where the coupled coefficients ^ (N a Nfo) are defined by 

1 2 


C kik 2 < N a N b> " C k 1 ( N a) C k 2 (N b><il m a’i2 m bl JM ab> 


+ '(- 1 > T c k 1 ' N b) c k 2 ( N a ) <Sl”b' J2 m al JM ab> 

with the proviso that only the first term shall be used when k^ = kg and 
T = 0. The coupled orbital potential energy matrix elements may now be 
expressed quite simply as 
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N'N t |V J |N a N b ) T = ^ 1 D k3k 4 n N b) C k[k 2 < N a N b){ k 3 k 4 IIV ' lk l k 2/ JT . 


k 2^1 k r*3 


where /kgk^HVllk^kgj are the conventional two-body coupled matrix 
\ / JT 


elements. 


NEUTRON-PROTON SYMMETRY 


For nuclei with N = Z, it is common to assume that $ has n-p 
symmetry; that is, for each neutron orbital n )' t ) in $ there is a 
corresponding proton orbital n_> -t_) identical except for isotopic spin 

a a . 

(and similarly for $’)• The functions n_) do not depend on r in this 

a a 

case, and hence the kinetic energy integrand may be written 


A/ 2 


t(9) = 2 Z I d mlm 

"i"a J 


V K\ p 


(1 >k 


>• 


Now, in general, when it occurs in v pn (0) the two-nucleon density 
matrix has only one term, 

< N k N bl p(2> l N a N b> = < N klP (1) l N a > < N bl p(1> l N b >. 

whereas in v (0) and v (0) both of its terms contribute. By rearranging 

P ** 

the terms in the series, however, it is possible in the case of n-p sym- 
metry to express the potential energy integrand in terms of the isospin- 
coupled density matrices 

(»in' b |p (2) |n a n b ) T = <n^|p (1 >|n a > <^|p (1 >|n b > 

+ (-l) T <n^|p (1 ^ n b ) <n{j|p^|n a >, 
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which are defined with the proviso that when n a = n^ and T = 0, only 
the first term is to be used. The final result is 


A/2 A/ 2 



R y (7r) SYMMETRY 

The operator which rotates the z-axis through an angle 6 about the 
y- axis is given by 


-iGJ_ 

R y (0) = e y. 

The special case 9 ■= tt is important, for the spherical basis functions 
transform in a particularly simple way when R y (7r) is applied: 

R y (7r)|km> = (-l)j~ m |k - m). 

In general, the orbitals In ) do not possess corresponding prop- 
erties . However, for nuclei with even Z and even N, it is often as- 
sumed that $ itself has R y (7r) symmetry; that is, for each orbital 
In > in $ there is a corresponding orbital |n_) which satisfies 

|N a > =V )|N a > l N a> = - R y ( ’ )|5 a>' 


As a result, the expansion of $ may be written 
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where, for convenience, we specify the subgroup {l N a >} of orbitals 
which serve to define $ 1 by the requirement that m >0 for each 
member |N a > . Clearly = Ry00$2> and one may also show that 
$1 = Ry(Tr)$2 and ^2 = ^y( 7r )]~ ^ The "rotated" state function $ 
then becomes 

~ -i0J„ +i7rJ 

$(0) = + e y e y $2 

= ^(0) + ^(0 - tt), 

arid, from the symmetry properties of the reduced rotation matrices, 

(jm|$( 0 )) = (jm|& 2(0 )} + (-i) m " k (jm|$ 2^ - 0)'). 

Consequently the integrand H(0) may be written 

h(0) = (:$’|h| 52(^)) + (-D K, “ K ($ , |H|i 1 (7r - e )), 

and, if <3? T also possesses R y (7r) symmetry, 

H(0) = 2 [($||h|$2(0);} + (-1) K, “ K ^^|h|$ 2(77 - 0))j. 



14 


In either case H(0) is of the form 11^(0) + (-1) K ~ K - e), and by a 

change of variable in the second term one may obtain 

Hj = r H 1 (0)[d^^(0) + (-1) J " K d£ t K ( 0 )] sin 9 de - 

2 Jo 

Since the quantum numbers m_ from In ) and |N ) will cancel in 

a a a 

pairs to give K = 0 (regardless of the symmetry assumed for $’), this 
reduces to 



Hi(0)djjr f Q(0)sin $ de, 


J = 0,2,4, . . . 

For the evaluation of V(0) it is convenient to employ a different 
specification of the manner in which $ is split into two functions 
$1 = R y (7r)$2 and This is accomplished by considering 

the expansion 


$ = (1 + 


A/2 

[l +Ry(ir)](2A(A - 1)]' 1/2 Y 


N„N. 
a b 


l N a N b>% N, 
a b 


and defining the subgroup {|N a N^>} by the requirement that M a b > 0 or 
M a b = 0 and m a >0. In this way only the rotation matrices d^j M ,(0) for 
M, M’ ^ 0 need be calculated. 

MISCELLANEOUS DETAILS 


By far the most time-consuming part of the calculation is the con- 

yrji 

struction of first the coupled coefficients ^ (NgN^) and then the coupled 

1 2 


matrix elements 


N’N’ 
a b 


V. 


N a N b 


This task is further complicated by 
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the necessity of keeping track of all the indices in an efficient manner. 

Probably the best approach is to order the orbital pairs (N a N^) first 
by isospin, then by parity, and last by M^. (Note that for a given J, 
only those pairs with |M a ^| < J will contribute. ) The basis pairs may 
be ordered in a somewhat similar fashion: first by parity, and then by 
J (such that I j i “ 3 2 1 — J + 32^’ T ^ us the P airs are arranged in the 
order in which they will be needed, and can be located by a single in- 
crementing index. 

Much time can be saved if the Clebsch-Gordan coefficients are cal- 
culated in advance, stored, and retrieved as needed. Ordinarily this 
requires some testing of the angular momentum quantum numbers to 
avoid dealing with those coefficients which automatically vanish. This 
testing is time-consuming, however, and it may prove to be more effi- 
cient in the final analysis to dispense with it entirely and simply multiply 
by the zero coefficient. The storage and retrieval problem can also be 
simplified if zero coefficients do not have to be avoided. Whatever the 
choice, it is worthwhile to check that min(j^, jg) — min(|m a | , |m^|), as 
these quantities can be calculated in advance when the pairs are ordered 
and later referenced by the same incrementing indices. 

Once the coupled matrix elements have been obtained, the integrands 
t (&) and v(d) are computed for selected angles by summing over the 
initial-and final-state orbitals. In the usual case where $’ = $, one can 
substantially shorten the calculation by making the replacements 


I - I 2 I f < N a' N a> 


N’N N’<N 
a a a a 


2 2 -* 2 




N’N’ NR 
a b a b 


N a N b^a N b 
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where €(a, /3) = 1 - 8^. Any suitable algorithm may be used for the 

numerical integration, but of course Gaussian quadrature is an obvious 

choice. It should be noted, however, that A (0) is a very sharply peaked 

function, so that some care must be taken in order to insure sufficient 

20 

accuracy. In a typical calculation for Ne, for instance, A(0) falls 
from 1. 0 at 0=0° to 0.07 at 0 = 45° to 0. 001 at 0 = 65°. In con- 
trast, h(0) varies relatively slowly over the entire angular range, 
changing by less than 50 percent. Probably the most satisfactory method 
of dealing with such an integrand is to base the Gaussian quadrature on 
polynomials orthogonal with respect to a weight function similar to A(0), 
but this is an involved procedure. Fortunately the zeroes of the Legendre 
polynomials are concentrated at the ends of the angular range, and this 
feature may be further enhanced by an appropriate change of integration 
variable before the quadrature method is applied. 
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